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ABSTRACT 


An analysis based primarily on computations is 
conducted to determine errors in a bombing system 
Similar to that currently under test and development 


by the Marine Corps. 


There are many parameters 20 the systen 
determining point of impact. Included are point of 
release and velocity at that point, air density,the 
mass and drag coefficient of the bomb, its cross 
sectional area, andthe wind structure. Most values 
are not known with the desired accuracy but must be 


estimated, say, from radar measurements, for position. 


This study estimates as much as possible the error 
in predicting the point of impact as a consequence of 
such factors as random uncorrelated noise in the radar 
measurements, errors in estimating density, errors in 
the wind structure, and errors due to initial 


conditions. 
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I. INTRODUCTION 


This thesis developes and uses many of the mathematical 
techniques and relations that are needed to analyze the 
important errors in a bombing systen. The AN/TPQ-27, 
described in Ref. 1 and Ref. 3, is the model for tke 


analysis. 


Precision bombing 1S a crucial means of supporting 
Marine Corps operations. [It is a matter of life and deatn 
that it be as accurate and effective as possible. To ensure 
this the Marine Corps carries out extensive analysis and 


testing of its systems. 


The complete analysis consists of two basic. parts. 
First there is a mathematical modeling of the system; that 
is, equations must be formulated to describe its behavior 
and ainterreactions of itsS components. These must be based 
on the natural laws that govern the moticns and on the 
accumulated experiences (some lore perhaps) with similar 
systems, The equations must be Simple enough to allow the 
necessary solutions and computations, but they must also 
Maintain the needed accuracy. This thesis takes up this 
part. Second, there must be extensive testing of the 
components and of the integrated system under conditions 
which Simulate the field as closely as possible, but with 
extra measurements made to isolate and determine the sources 


of Significant causes of error. 


In this study the author aas set up the basic 
Mathematics for error analysis and carried out computations 


to show the typical effects of errors caused or observed at 


ae 








various stages. 


Bees convenient to consider three classes of errors, 
depending on the time they are observed. Different analysis 


is required for each type. 


First are the errors which are due to inaccuracies in 
the equations of motion during the period of fall. These 
may come from sources such as’ errors in estimating the drag 
coefficient, and variations in the air density and the wind 


structure from the assumed values. 


Second, there may be significant errors in the position 
and velocity at the time of release. These may come fron 
sources such as wind gusts near the time of release, faulty 
plane response, errors in the filtered estimates of position 


and velocity, and release or separation problens. 


Third, there may be errors in the radar measurements and 
transients introduced in various ways. The errors may be 
inherent in the radar, they may be due to the changing 
aspect of the plane, to maneuvers, and/or to the filtering 
equations whose function is to reduce the effects of random 
errors. We will Know only statistical properties of some of 
these, and for these we can only determine associated 


statistical properties of the associated miss distances. 


A number of estimates were made in two or more ways, for 
comparison of the methods and to act as ae check. For 
example the effects of variations in initial velocity were 
estimated in three ways: (1) by recalculating the 
trajectory, (2) by using variational equations, and (3) by 
using the adjoint or costate equations. For small changes 
the latter two are more accurate; for large changes, the 
former is better, but all are adequate over a considerable 


range. Estimates of errors in the filtered values of 





position and velocity were also obtained in two basic ways: 
(1) by an analysis of the solutions to the difference 
equations that are used for filtering, and (2) by a 


numerical simulation of the solution, on the computer. 


We need to make some remarkS avout notation and 
coordinate sets at this time. We will use an x-y-z 
Cartesian coordinate set fixed with respect to the earth. 
The program haS a Subroutine which corrects automatically 
for the earth's rotation so we can consider the set to be 
Newtonian. The 7z direction is the vertical. The plane is 
assumed to be flying horizontally, and its velocity (with 


respect to the ground) defines the y-axis. 


Derivatives with, respect to time are indicated by a dot 
over a variable, &¥=s hy [at , etc. In the text, however, 
this is denoted by xdot, and a similar notation is used for 


other variables. 


The ratio "po is used in the equations as the ratio of 
the airc density at altitude z to the air density at ground 
level (We use an exponential value for this.). In the text 
thes= terms are written rho and rho(0) respectively. In the 
Same Way variations such as SX are written as DELx in the 
text. 


It was not feasible to write subscripts, so that 
subscripted variables such as the drag coefficient, the x 
component of the wind velocity, the y component of the wind 
velocity, and others, were written Cd, Vwx, Vwy, etc., in 
the text; it is felt that this will not be confuSing. 


The notation was chosen to be in conformity with that in 
Ref. 1 and Ref. 3 as much as possible. Also, the terms 
impact and fall are both used interchangeably to denote 


conditions at the instant the bomb strikes the ground. 





II. ERRORS DUE TO CHANGES IN DENSITY, WIND STRUCTURE, AND 
VALUE OF DRAG COEFFICIENT 
When tne ordnance is released it follows a path 


determined by natural laws. We approximate these as best we 
can to predict the point of impact. However, the equations 
of motion contain certain parameters and functions which may 
not pe known, and the planned trajectory has a point of 


impact different from that calculated, and a miss results. 


The effects of two of these parameters are investigated 
here. The first of these is the wind. Let us assume that 
we have a gcod estimate of the wind at the pcint of release 
(called the wind at altitude), because we know the velocity 
of the plane (from the radar), its airspeed, and its 
heading. However, the wind varies with altitude, and the 
structure of this function may be unknown under combat . 
conditions. When testing, we can approximate this function 
by making extra meteorological measurements. Then by using 
Various approximations to the wind, based on the wind at 
altitude, we can calculate the associated miss distance. An 
approximation that is commonly used is that the wind is 
constant or uniform, that constant being some fraction of 
the wind at altitude. The lore indicates that eight-tenths 
1s a good fraction; higher values, approaching one, have 
been suggested [Ref. 4]. Under the conditions of our 
example eight-tenths seems like a good choice. However, the 
component of the wind normal to the wind at altitude often 
causes a iarge error that there is no way to correct knowing 
just the wind at altitude. Only a further study of wind 
structures will show when bombing is ineffectual without 


more detailed knowledge of the wind. 
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A second factor whose effects we can estimate readily is 
the drag parameter. There is a single drag parameter, which 
involves the air density at some reference level, the cross 
sectional area of the bomb and its mass, and a drag 
coefficient, Cd, depending primarily on its Shape. We can 
estimate the effects of errors in this parameter directly. 


There are many other sources of error for which no 
Simple analysis exists; some of thes? are discussed later in 
this chapter. 


A. EQUATIONS OF MOTION FOR A PROJECTILE 


Referenc2 1 describes in detail the equations involved 
in the ballistic integration. These are th2 classical 


equations of motion for a non-lifting projectile: 


Ke Alm p GV (KY) 
We (/6 ) 
G = -Alm p C, CY (g-V y) 
z= -9-Alm p CV 2 We) 


( 


Z 


where A is the cross sectional area of the projectile, mis 


the mass of the projectile, gis the acceleration due to 
gravity; rho, Cd, xdot, ydot, zdot, Vwx, and Vwy have been 
previously defined; and V, the airspeed, is: 


- V(x - . ,- *e° (a) 
VVRLTG YEE 


Equations (1a), CTD, and (lc) were integrated 


humerically to obtain final’x, y, and 2Z coordinates and 
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hence, to estimate the bomb impact point. 


B. BASIC ASSUMPTIONS AND CONDITIONS 


It was decided to use values for various parameters that 
were selected arbitrarily in a reasonable range, rather than 


treat a specific missle, whose details would be classified. 


i. Altitude and Airspeed 


= 


In all cases the aircraft was assumed to be flying 
straight and level at an altitude of 20,000 feet with an 
airspeed of 350 feet per second. In those cases where a 
wind structure was considered tne initial velocity (with 
respect to the garth) was altered to take into account the 


wind velocity. 


2. Drag Coefficient and Constant Values 


aS ee SS ye ee SS er SP SSS a Pee ep ee 


The drag coefficient, Cd; the cross sectional area 
of the projectile, A; the mass of tne projectile, m; and the 
air density at ground level, rho(0); were assumed for the 
Majority of the work to take on values such that the 
important ratio C=(A/m) (Cd) (rho(0)) had the value 0.000025. 
This figure was considered to be a good approximation to 
some typical ordnance tyves without reguiring the thesis to 


be classified. 


3. AL 


ir 
fw) 


easity 


It was assumed the decrease in air density with 
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altitude was given by the factor exp(-z 1n2/18000). This 
function was compared with the ratio given in Ref. 1, and it 
was felt that the exponential expression was satisfactory 


(and simplified programming). 


Reference 4 details wind structure data obtained by 
Camp Pendleton Marine Corps Base in 1975. The data was 
collected twice daily on 25 March, 28 March, and 12 April. 
Wind velocity in knots and wind angle were determined at 
2,000 foot intervals from ground level to 36,000 feet. All 
reference within this work to a "variable wind structure" 
applies to the data in Ref. 4 collected 25 March, 1975, at 


11:00 unless otherwise noted. 
CA COMPARISON OF SEVERAL WIND CONDITIONS 


Several trajectories were computed based on two actual 
Wind structures. In each case directions were chosen for 
the plane corresponding to a headwind, a crosswind, and a 
tailwind. In each case trajectories were computed with a 
uniform wind equal to the wind at altitude, and one equal to 
eilght-tenths wind at altitude. In two casesS various other 


unirorm winds were also assumed. 


In all calculated trajectories the following equations 


and terms were used: 


x = CL (32) 


(34) 


——~< 0 
TT 
~ 
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All trajectories were calculated uSing a modified Euler 
Method [Ref.2] with a time step, DIT=1.0. 


Since the wind velocity components in Ref. 4 are 
provided at 2,000 feet intervals only, a linear 
interpolation scheme was devised to compute wind velocity 


components as the bomb fell. 


Three trajectories using the wind structure at 11:00, 25 
March, 1975, were then computed for three conditions: a 
plane flying into a headwind, a plane flying with a 
crosswind, and a plane flying with a tailwind. Next, wind 
at altitude was considered as the constant wind velocity 
throughout aS opposed to using a variable wind structure, 
and the three trajectories computed again. Finally, the 
three trajectories were computed again asSuming the wind to 


be uniformly eight-tenths that at altitude. 


the trajectories based on the variable wind structure 
were used as trajectories of reference and the errors due to 
the approximations calculated. Corresponding trajectories 
using the stratified wind data of 25 March, 1975, 14:00, 
contained in Ret. 4 were then calculated under the 


conditions described above, for a comparison. 


Tf the wind at altitude is known, and is all:that is 
Known about the wind, the only uniform wind we may 
reasonably choose is one which is some fraction of the wind 
at altitud?. It was felt that the change in impact point 
due to changing this fraction would be approximately linear. 
To check this two trajectories were calculated using 
nine~tenths and zero respectively, as fractions. The 
resulting trajectory uSing nine-tenths led to a point of 
impact which was approximately the midpcint between the 


points due tc fractions of @ight-tenths and cone, indicating 
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the change waS approximately linear. A typical plot (Figure 
1) Showed that in the example considered, the value of 
eight~tenths waS near the optimum. It is felt no other 
value would result in a significantly smaller error in any 


of the cases considered. 


While Ref. 4 points out that fractions of wind at 
altitude closer to unity are better for certain bomb loads, 
the results indicate that for these wind structures, 
eight-tenths is near the best choice for use in the absence 
of stratified wind data. The resulting error was near the 


Minimum in all cases for this choice. 


For these wind structures there was always a sizeable 
error when a constant wind was assumed, due to changes in 
the direction of the wind with altitude. More analysis with 
different wind observations is needed to establish the 
validity of such a universal factor and the general errors 


due to the varying direction of the wind. 


Some additional results of the analysis are contained in 
TABLES I and II. 


In the next chapter another method is given for 
estimating the effects of the changes discussed above, using 
methods like those developed by G. A. Bliss for corrections 


in artillery. 
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Figure 1 
A typical Plot of the Errors From the Variable 


Trajectory Impact Point (in Ft.) When Various 
Fractions of.Wind at Altitude Are Assumed 


Nd 





Qe De OWN ATA FROM 25 MARCH, 75, 11:00 


_ type wind vel x y x error y error 
wind used 
ened variable 67.754 15519.350 N/A N/A 
Gay. wind at -.028 15629.180 =O1.6 6, 
alt 
ne weal ~8 wind -.023 15499 6.550 -67.8 -20 
at alt 
eross variable -409.653 10994.860 N/A N/A 
eross wind at = ee oe 1090222150 =—102,9 -{2.1 
alt 
cross 73 wind -412.139 10929.470 -2.5 -65.4 
ab ale 
head variable -65.814 7488.727 N/A N/A 
head Wane ea t O20 Te ou 65.5 ee 
au 
head .8 wind O70 15072962 6523 1949 
at alt 
Table I 


X and Y Errors (in Ft.) for Wind at Altitude 
and .8 Wind at Altitude Trajectories Versus 
Variable Trajectory Under Various Wind Conditions 
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type 
wind 


tail 


tail 


tail 


cross 


CcrOos 


Cross 


head 


head 


head 


wind vel 
used 


variable 


wind at 
alt 


~8 Wind 
at alt 


variable 


wind at 
alt 


.8 Wind 
a alt 


variable 


wind at 
ale 


-on vine 
at alt 


STRATIFI“£D WIND DATA 


ROM 2 DO MARGH st. 


x 


57.813 


-.028 


-.023 


=-422,153 


-566.438 


-449.604 


-56.073 


0.0 


Oe 


y/ 


1535 5.250 


16010 ..070 


15868. 360 


10874 .430 


10806.460 


TOS to. 00 


1142.984 


6995.461 


7139.324 


Mev 1WIE 


x eri Or 


N/A 


-57.8 


-57.8 


N/a 


-144.3 


-2/.5 


N/a 


Dono 


a0) 5 © 


14200 


y error 


N/A 


(AS 


N/a 


-68.0 


-59.0 


N/A 


- 148.0 


-3./ 


X and Y Errors (in Ft.) for Wind at Altitude 


and 
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.8 Wind at Altitude Trajectories Versus 
Variable Trajectory Under Various Wind Conditions 





D. COMPARISON OF VALUES OF THE CONSTANT C 


The value C= (A/m) (Cd) (rho(0))=0.000025 was, as stated 
earlier, used for the majority of the computation. [In this 


section this parameter was varied to determine the effects. 


Trajectcries were computed with various values of C, 
namely OO, .0000125, .000025, and .000050, to determine the 


associated change in point of impact. 
As indicated in TABLE III, the function is nearly linear 


over a wide range, and hence changes in point of impact 


could be estimated by linear approximations. 
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sV0O025 


200050 


OOOO 125 


Iigo2. 190 


OO TOO 


I SV2 ee Ree Sl, 


[AS ara oe) 


abi ex en 


CHANGE IN 
IMPACT POINT 


N/A 


194.97 


425.76 


Si On. 


Change in Point of Impact (in Ft.) as Result 
of Varying the Constant C (Trajectory with 
C = 0.000025 is Nominal Trajectory) 
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From the point of view of control the time of release of 
the bomb is the crucial instant. Conditions at this time 
determine success or failure, and nothing rurther can be 


done to improve bombing accuracy. 


In this chapter a study is made of the changes, or 
errors in the point of fall, as a consequence of changéS OF 
errors in initial conditions. By initial conditions are 
meant the conditions at the instant ot release from the 


plane. 


The effects of changes in the initial values of x and _ y 
on the point of impact are obvious. The effects of changes 
in u, V, Ww, and z are not so cbvious; these effects are 
estimated in three ways and the results compared. The first 
way 1S to take a typical trajectory and then recompute the 
trajectory several times, | each time changing one initial 
condition by a "small amount". The second way is to use the 
Variational equations. These were integrated, uSing the 
appropriate initial conditions in order to derive the 
differential estimates of the changes in the point of 
impact. The third way is to define the adjoint or costate 
equations. These are integrated backwards, with appropriate 
terminal conditions, to define the so-called influence 
coefficients, which may .be interpreted as the partial 
derivatives of the final conditions with respect to the 


initial conditions. 


Compariscns of the results were carried out with three 


reference or nominal trajectories. In one trajectory there 


PaP 





was asSumed to be no wind; in the second the wind was 
assumed to be uniformly eight-tenths that at altitude; in 
the third the wind function was assumed to be that observed 


at Camp Pendleton on 25 March, 1975, at 11:00. 


In all cases all of the methods gave vary similar 


results. 
A. VARIATIONAL EQUATIONS 


The variational eguations are derived as follows. 


Consider the eguations of motion: 


az -CeVCu-Y ) (Ja) 


L/X 
oO 
v- -CeV(-Y. ) (1) 
st 

Pa) 

w= -CpVw - ¢ (Ic) 
where u, V, W, C, rho, V,Vwx, and Vwy have been previously 
defined. Then let us consider the first-order effects if we 
change the variables by small amounts in some way, but 
require that the new variables satisfy the eguaticns of 


motion. The first-order terms in the change, DELx, etc., 


must then satisfy the equations: 


Sx = Su (2a) 


Sy=SV (24) 
£2-Sw Coc) 
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These are called the variational equations. It may be 
noted that they are linear but with variable coefficients 
whose values must be obtained from the original trajectory. 
Jf the variations in the wind are zero, the equations are 


also homogeneous. 


1. Effects of Changes in vy 


ep te 


As a first example let us consider the effects of 
changes in the initial speed, v(Q). Let us take as a 
standard or nominal path one with no wind. The values of Z, 
u, V, W when they occur in the variational equations are 


from this path. 


Let us consider DELv(0)=1, and all other variations 
to be zero initially. We will integrate to the terminal 
time tf, the time of fall on the original fath. We find 
that DELy(tf)=31.30. However, DELZ(tf)=1.34, and hence the 
differential change in terminal time (denoted Dtf) is 
Dtr=-DELz (tf) /w(tf)=0.001. The dirferential change in y 
(denoted Dyf) then requires a secondary term due to the 
change in terminal time: Dyf=DELy (tif) +v(tf) (Dtf)=31.63. We 


see that the change in terminal time has a secondary effect. 


The variation in x in this case is zero. 


2. Effects of Other Changes 


In a Similar manner the effects of changes in other 
1hRitial conditions were studied. The variations DEILu(0), 
DELZ(0), and DELw(0) were each Separately set equal to one 
initially, while all of the other variations were set to 


zero initially; we integrated to terminal time, tf, and 





computed Dxf and Dyf in each case. We then computed the 
error in impact point due to the differential changes in x 


and y. 


Identical calculations were performed considering 
the three separate nominal trajectories, and as indicated in 
TABLE IV, differences in figures among the three 


trajectories are quite minimal. 
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NOMINAL #RRORKR WHEN #RROK WHEN #RROR WHeAN ERROR WHEN 
RAJ. DELV(0)= v(0) DubZ(o)= Z(0 ) 
eo INCREASED rae ENCREAS EE 
ono) ORS 
no wind S\eG95 515.958 528 bar) 
~8 wind Slr ed SO eG Ale AS Ip 
eV voees Wear 
variable Seo. 4 Sen Or! Paes ores 
wind : | 7 
HRROR WHEN s&RROR WHEN cRROR WHEN S#RROR WHeN 
DELW(0)= W(0) DELU(0)= U(0) 
ee INCREASED 1.0 TNCREASED 
Vos) (ropes) 
no wind O75 55 96. 505 32.867 328.646 
~5 wind Soe SA oe B2.109 623.0 
eae 2G 
variable S506 85.60 D2 elco D2 a4 
wind 
Table IV 


Comparison of Errors in Impact Point 


(in Ft.) When Errors in Various 
Initial Conditions are Assumed 
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B. THE ADJOINT SYSTEH 


The adjoint system For the equations of motion was 
derived, and the influence coefficients analyzed aS a means 
of cross checking the results of the variational equation 
analysis. The nominal trajectory for this ccmparison had no 


wina. 


1. Derivation of Influence Coefficients 


— —_— a a ee ee ee ee — ar 


Let us introduce six new unspecified variables 
P1,P2,P?3,P4,P5,and B6 to be defined later. The hamiltonian 


for the system is defined: 


= Put hye hw-pCV | (3) 


LAGU) eG, )eeul- 2 y 


wher eo Vy WW, rho, C, Vwx, and Vwy have been previously 
defined. Now, let us always choose the P's to be solutions 


to the differential equations: 


re) 
Pe, % dH (4) 
C 

dX, 
where xX1=x, xX2=y, xX3=z, x4=u, x5=v, and x6=w. 


Then 


-_ -dtH (Sa 
J = AX SOE ) 
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le 


a 0 uo 


ROO 


f- 8 


dy 


2 fee - Fst cv(2le-, )e“(- U,, Je eu) 


3 cla l- Ve.) lM, + eer J 
= -_ - rep QiVh + ole-bo)+66-Y, eeu] 


oe 
V 


“hep CVe, wader 





seks LeBel) 


Hap Clverlak-l, va Ke bes] 


(} 


= teepvbep Phe) GU, Jot, tel 
8 Snpcltay deel, bed ff 
= -f+pCVE+ ee lel-Y, je bl-Y buf . 


(54) 


(5<) 


(5d) 


(Se) 


(SF) 


(59) 


(54) 


(5 /) 


(Sj) 








This does not define a solution yet. If we define 
values for all six components at some particular time, a 
particular solution will be determined completely, though we 


may have to integrate numerically to obtain it. 


The adjoint variables and the variations satisfy an 


‘important identity 


6 


a . 
t=/ . 
which iS easily verified, whenever the DELxi are solutions 

to the homogeneous variational equations. 
This relation is the basis for obtaining what are 
sometimes called influence coefficients cr sensitivity 


factors. 


To illustrate how these are obtained, let us derive 
the senSitivity factors for the final value of y, cr x2. 
First, let us observe that if we integrate equation (6) fron 
zero to ti we get 
: (7) 
> Pa = (- (2 SFX, 

7 2 A c 0 


: é 
Ui 
Now, let us choose a solution to the adjoint such that at 
time tf the coefficient of DELy=DELx2 is P2(tf)=1, but all 
Other Pi(tf) are zero. To generate this solution we must 
integrate backwards to t=0. For this particular solution to 


the adjoint system we get then from equation (7) 


6 
(4) =(& oe " 


t=O P 
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That is, for this particular solution the components Pi(Q) 
yield the influence of DELxi(0) on DELy(tt), and hence the 


hame. 


Several comments may be made. There is one such 
solution vector Pi for each variable x, y, Zz, uu, v, and w. 
These were calculated for their corresponding variables. 
The values obtained in this way for the final variations 
should be the same aS those obtained from the variational 
equations except for errors in numerical integration due to 


EQundofE, truncation, etc. 


2. Analysis of Results 


Now, the influence coefticlents are actually the 
various partial derivatives of the argument at final time 
with respect to the arguments x, Y, 2, wu, Vv, and w at 
initial time. For instance, the influence coefficient 
associated with x defines a six component vector whose 
components are the partial of x at final time with respect 
to x at initial time, the partial of x at final time with 
respect to y at initial time, the partial of x at final time 
with respect to z at initial time, the partial of x at final 
time with respect to u at initial time, the partial of x at 
final time with respect to v at initial time, and the 
partial of x at final time with respect to wat initiai 
time. The values of all these partial derivatives or 


influence coefficients are given in TABLE V. 


These coefficients serve as a cross check for the 
analysis of the variational equations. For instance, the 
Variational eguations with DELu(0)=1.0 yielded a value of 
DELx at tf of 32.867 while the adjoint indicates that the 
partial of x at final time with respect to u at initial time 


is 32.870, which agree to four figures. Likewise, the 
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variational equations with DELv(0)=1.0 yielded a DELy at tf 
of 31.297 while the adjoint indicates that the partial of y 
at final time with respect to v at initial time is 31.298. 


Other values were checked Simiiarly. 


In conclusion, we found that the three mnethods of 
calculating the effects of initial conditicns on the point 
of fall: (1) recomputing the trajectory, (2) uSing the 
variational equations, and (3) using the adjoint are 
satisfactory within the range of interest. We did find that 
the differential method augmented by the correction fcr the 
change in terminal time of fall is slightly more accurate, 
but the increase is really not Significant. Finally, we see 
that in a computer age such as today the effects can just as 
easily be calculated by changing the initial condition and 
programming the trajectory again, although the use of tae 
variational equations iS still an important tool for 


analyzing such problems. 
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Xe /al New Y0;%e, LoVe, “W/o » le Ja Vie! Yo, To, ke, Vo,Wo 


lee) O.0 
0.0 en) 

O 0 702 
BZ or {He 0.0 
0.0 Ble Sie! 
O20 22D 


Oz¢ [dC ) Iyer Yo 2o,Ule Vo, th dug [oC Vee YU, Zo, Uo, Vo, Wo 


ae 0.0 
hoe 0.0 
920 0.0 
0.0 crow 
1.338 a8 
28.983 ee 


Aale ) bee, yo 20, tee, Vo, te We /d¢ Vicar us, Ewe 


O20 OO 
On 9 20 
O004 SOLON 
O20 O. 8 
.616 064 
~147 O07 
Table V 


Adjoint Sensitivity Factors 
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C. COMMENTS 


The problems of Chapter II can also be solved by the 
methods of this chapter. Indeed, Bliss developed the 
techniques associated with the adjoint to handle just such 
problems, while he was at Aberdeen Proving Grounds during 
World War I [Ref. 7]. 


Let us consider the effects of varying the wind profile 
from its values on some nominal path. The variational 
equations, (2a) through (2f£), apply directly. The values of 
Z, U, Vwx, etc., are those from the given nominal path, and 
the values of DELVwx, DELVwy are the values on the new path 


Minus those cn the nominal path. 


If the adjoint equations are used, the equation for, say 
DELy (tf) , reduces to 


Sy4)= ‘C/ Venp C thm /i2000) (9) 
[Be Vino, 4 BOY ek, Jee 


where the solutions to the adjoint are those associated with 
y, for which P2(tf)=1, and the final value of the other Pi 
are zero. The evaluation of the integral is simplier than 


integrating a system of differential equations. 
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IV. ANALYSIS OF FACTORS INFLUENCING AIRCRAFT POSITION AND 
VELOCITY ESTIMATES 


The effects of some errors in radar measurements of 


position on the filtered values for position and velocity 


are investigated here. Also included in the study are 
errors due to some simple Maneuvers and associated 
transients. Reference 1 provides the equations used in tne 


AN/TPQ~27 system to predict and filter the aircraft position 
and velocity. The analysis is generally carried out for a 
Single coordinate, y; by definition xdot is zero and 


ordinarily zdot will also be Zero. 
A. RANDOM UNCORRELATED NOISE 


The first error analyzed is that due to random 
uncorrelated errors or noise in the radar measurements. The 
variance of the resulting errors was determined in terms of 
the variance of the noise, in two ways. The first was 
purely computational. A set of random numbers, v(k), with 
mean value zero, variance one, and gausSian distribution was 
generated. These were used like measurements Or 
observations in the filtering equations. Sequences of 
filtered estimates for y and ydot were generated and their 
means (Which were near zero) and variances were calculated 
On the computer. The second meéetnoad was to analyze the 
filtering equations to estimate the variances of the 
filtered estimates of y and ydot. 
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A routine for generating random numbers was called 
from the computer library. It generated random numbers 
V(K), k=1,2,....etc. These values were used as a forcing 
function, y(kK), in the filtering equations. The position 
and velocity, y and ydot, were initially set to zero, and 
then their filtered values calculated iteratively. 
Numerical estimates were then made for the standard 
deviations of the filtered values of y and ydot using those 
calculated values as the sample [Ref. 5]. It was then found 
with this strictly numerical analysis that the standard 
deviation of the sample of filtered values of y was 
approximately .245 % and the standard deviation cf the 


Vv 
sample of filtered values of ydot was approximately -1187°. 


- _ y wa y AT (/) 
m, m=! m-tm-l Mm-l, m=! 

a A 

a a (2) 
mMm-| m-l, m-! 4 


The filtering equations are: 


ye 7 rey. g ) G3) 


my m-t 

A A 

o oO A (4) 
: + # ( - ) 

yey (nS Mae, 

mm mm-| ™ ha FE fa | 

where "7 is the estimate of y(tn) nade at time tm; the 

corresponding estimate for ydot is y ; alpha=.1, 


beta=,0052, and DT=,125 [Ref. 1]. Now, the “Substitution of 
the right side of (1) into (3), and the substitution of the 
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Light sides of (1) and (2) into (4) yieids: 


at +a (yf -4 ae} (5) 


A 
+ o 
/, : / ~f ~/ 
,™” M-l,m=} ml, mel m-l,m~] 4m! 


A-!, m=} 


(I-«) 4 + a7(1-«) y + ty 


+h -£ og “) ai 


3 


3 ™Mos 
\ 
yNos 


These are the eguations for the filtered estimates of y and 


ydot, in terms of the observations. 


Next, the following substitutions were made: 


y =-4 (7) 
mm ” 
; 7 t (g) 
mm ‘at 
A 
or ae (9) 
",™ 


then (5) and (6) become 


Dani)? +d tay (0) 


m~/ 


F. 7 + (-#) a ud 4 | “ 
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One of the variables, say, zeta may be eliminated as 


tollows. Shift the indices in (10) and (11) by one, to get 


1 > U-«)7 + /-«) P ea (/2) 


y 2 -# 77 0-8) ¥ i el (13) 


When zeta is eliminated among these equations the equation 


“4 = Q-«-4)/7 -¢-Y7 lad eee (/#) 


If we set 


A = a-X-f (/5) 
B= - (i-«) z 


C = £-X (17) 


this becomes 
VO - AH +BY +aeVu CV (/8) 
4d m~/ m-2 Sh: i ta 4 


in the same way the equation for zeta may be 


oktained, 


eral wre AG + #V + OY (9) 
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where 


D=-6 . (29) 


In order to estimate the variance of eta we need to 
express eta in terms of the uncorrelated observations. To 


do this let us shift indices in equation (18); 


1 2 es a ‘or 6 CY (2/) 


™~/ 


and eliminate nak, from (18) with this; 


i -(A+b)0 + ABY ta (22) 
+ (fat C)V va ACV 


Repeated substitutions of this type lead to a series which 


may be written 


a 


™. 


«V+ Gar CW +a tr acy. _ 
+ (he +A ap A ee eee 
< m-3 


lV = 


kay 
m 
SS 
+ 
i. 
~‘. 
mh 


om ome m-2 
es, EE es 
et eee 
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The coefficients A, B, and E are defined by the 


recursion formulas and initial values indicated below. 


“= AR t @ 


an -/ m-/ 


(25) 


6 . BA (2é) 


ama 


E = H x Ft A & (27) 
(2€) 
O C27) 


E=o (30) 
O 


If we sguare both sides and consider the expected 


values we get 


2 i z (3/) 
Sz - C- 
(fe). 


A Similar operation for zeta leads to a recursion 


formula 


pae Af +f O 
44d pry am-/ 
where the A's and 8's have been previously defined and 


a : (33) 


We then ger 


g (ZF ‘)o os 


(32) 
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The series represented in equations (31) and (34) 
were summed iteratively. Calculated standard deviations for 
the sample of filtered values of y and of ydot were found to 
be .275 ve and .095 se respectively. The two methcds of 
analysis yielded values differing by approximately eleven 
(11) per cent for the standard deviaticn of y and by 
approximately twenty (20) per cent for the standard 
deviation of ydot. The difference of twenty (20) per cent 


ils somewhat disappointing. 
Be. PERLODIC ERRORS 


Measurements are also often sunoject to some scrt of 
periodic errcr, or noise. This was approximated by a sin 
wave and tne effects on the filtered values of y and ydot 
estimated. Various freguencies were considered and the 
resulting error in position and velocity was determined for 


each. 


The filtering equations (10) and (11) were of the 


Porn: 


1? 247 wa a? eA C35) 
~ ay mn 
y 


ax) 


© “87 ir cf + oy (326) 
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where 


B« (-4 Gy 


AS IS ‘ (38) 


The errors in measurements were expressed in complex 


form 
ac 
m 
where theta is the frequency. Multiple values of pi, 


differing by factors of two, were used for theta, and the 


amplitudes of the resulting steady state solution were 


calculated. 


The solutions were assumed to have the same forn, 


MOC (40) 


(ia 
£-Ge"™ : (1) 


and these were’ substituted into equations (35) and (36). It 


is convenient to shift the index by one in (35) and (36). 


YA - aff raP anal | (42) 


mt / 
P Z -8/ t ot +OY (43) 
WN |} ~ ty Ott ; 
We aiso shift the index in (40) and (41) by one. 
(nm HJ) Oc y 
A ae A (4) 
m t/ : 
Gg 
Or > ae (5) 
m t/ é 
Substituting (44) and (45) into (42) and (43) respectively 
yields 
C2 = ac #20, + oe? (46) 
Bé ; 
Cie “2-86, +4, +fe°* , (7?) 
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Now, (46) can be written 


(e*2)G- 2G = xe (48) 


and (47) can be written 


SC, + ‘ce +2 


Equations (48) and (49) are two equations in two 


be” (49) 


unknowns, C1 and C2. The solutions are 


C= (Fe BU 2 $é) (so) 


where 

Fe X Cos 26-244 Cos P@taroos 8 (51) 
Ge Si7 28-428 SI2n BAP Sin O (52) 
Hz 2ae+aPt cos L8- (are)tos 8 (53) 
S= SIN R8- 2517 8-1 5/7 8 C54) 


an gb #5¢) (55) 


where 
P= Pcos 28- 2Ff cos9- Pax Cos G (56) 
Pe S17 19-ak sing ~ PX Sing ; (5?) 
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Then the amplitudes of the Sin waves are 


[ ¢, | = [rueyys? — 
Ic. | - |/ P+ PIM 4S | (59) 


The magnitudes of C1 and C2 were calculated on the 
computer for seven different periods, and the results are 
shown in TABLE VI. 


These results show that when the errors in the 
Signal are Sinusoidal, the errors in the filtered values are 
also sinusoidal. The amplitudes of errors in the filtered 
values are much Smaller than the measured values when the 
period is small. For large periods the amplitudes of tne 
errors in the filtered values of y are appreximately equal 
to the measured values of y. Periods of order eight (8) to 
sixteen (16) seconds appear to be the critical periods, 
where maximum errors in the filtered values of y and ydot 


occur. 


As an example of the significance of a Sinusoidal 
error in the measured value of ydot, suppose we have a 
Sinusoidal error with an amplitude of ten (10) feet per 
second and a period of eight (8) seconds. This would result 
an an error in the filtered value cf velocity of 
approximately four (4%) feet per second, which would cause a 


Miss distance of approximately 125 feet. 


AA 








ee PERIOD (secs) C, f ) | coat 


A 25 051 024 
7/2 ne HORS $052 
mf 4 1 ~ 154 056 
7/3 2 265 me 
w/16°- 4 5 Sey Reo 
mw /32 3 eOs4 . 392 
wr /64 16 1.239 Moke. 
17/128 32 109 lige 
Table VI 


Calculated Amplitudes for Various Periods 
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C. EFFECTS CF TURNS 


Two types of maneuvers are analyzed here. 


First, during the final stage the plane can be expected 
to make a seguence of corrections, Which we may expect to be 
Small, Sharp turns. Hence, we need to know the errors 
introduced in the velocity and position estimates by such 
Maneuvers. AS a typical turn we considered a plane flying 
at 350 feet per second turning at a rate of three (3) 


degrees per second for one second. This yields values 


7al* 7.15 /6$ t-* ) KEf (C04) 


"{ 


¢ Ck) + le3fe k 5 kS8 (08) 
y CK) 


cy) 


normal to the initial motion. The values along the initial 


QF + /¢.3 /¢ (k - @) Pie (60<) 


t\ 


\l 


/&. 3 Fe oe | (ood) 


course are not changed significantly. 


The values of y and ydot and the filtered values of y 
and ydot were calculated, and the errors checked. The 
Maximum errors were of the order .0005 feet for y and .0008 


feet per second for ydot, which are negligible. 


Second, we considered a plane in a constant turning 
movement. This is an unlikely maneuver for an aircraft to 
be performing, Since the airplane flies on a straight line 


except for minor adjustments in the precision guidance 
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Stage. For this reason higher errors were expected. 


We calculated the errors in the fiitered values of y and 
ydot, and the corresponding miss distances as a result of 
these maneuvers. It was found that errors in the filtered 
values of y caused a miss distance of approximately 100 
feet, while errors in the filtered values of ydot caused a 
miss distance of approximately 1100 feet. This is a 
Significant miss distance, but as stated earlier, it is not 
surprising since a constant turn by the plane should only 


occur in the preliminary stages. 
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Several conclusions and some comments are presented here 


as a means of summarizing this study. 


This thesis had several purposes, focusing on the desire 
of the author to study mathematical methods and techniques 
that have been and/or will be useful when applied to a 
problem of interest to the Marine Corps. The techniques and 
methods used here are all known, but there is no place where 
the combination of numerical mechods, ballistics, 
differential techniques, difference (filtering) equations, 
and statistics is available; they generally are scattered 
thrcughout various fields. 


A. CONCLUSICNS 


Most of the results obtained in this thesis were, in a 
sense, expected. The range of validity of linear 
approximations to changes was surprisingly large, however. 


In our very limited checking, the fraction 2ight~-tenths 
of wind at altitude, actually given to us as part of the 
"lore" was found to be near the best choice for use in the 
absence of stratified wind data, Its validity, however, can 
be checked only by some actual testing and studying of 


actual wind structures. 


we also found that errors in various initial conditions 


result in errors in the impact point of the bomb which are 
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linear over a fairly wide range. It does not seem to make 
much difference which technigue we use to estimate these. 
Under the conditions of our trajectories an error or change 
of one foot per second in horizontal velocity led to an 
error or change in point of impact of about 32 feet, roughly 
independent or the direction of the error. A value of azdot 
of one (1) foot per second at release led to an increase in 
range of about 2.26 feet. The miss distances were found to 
be nearly linear over a wide range of values for the initial 


errors. 


The analysis of the difference equations was used to 
show how a random error in measurement leads directly to a 
random error in the point of impact. From the standard 
deviation of the radar measurements we derived the standard 
deviation of the estimates for y and ydot and their 
correlation. With these and the Sensitivity coefficients we 
obtained the standard deviation of the associated misses. 
Under the conditions we found that if the radar has a random 
error in measurement with standard deviation io , then there 
Will be a resulting error in point of fall with standard 
deviation } = 3 - For example, if the radar noise has an 
ras (root-mean-sguare) value of ten (10) feet, the point of 
impact will have an associated rms value of thirty (30) 
feet. Most cf this results from the error in predicting 
ydot, though y and ydot have a fairly large positive 
correlation, and hence the error is larger than if they were 


independent. 


Also we found that the presence of a sinusoidal error in 
the measurements resulted in errors which are sin waves of 
varying amplitude in the position and velocity estimates. 
We found that a sinusoidal error with an amplitude of ten 
(10) feet and a period of eight (8) seconds will cause a 
periodic error in the point of fall with a maxinum amplitude 
cf 125 feet and an rms value of about 96 feet. The actual 
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error depends on the particular time of release. Sinusoidal 
errors with periods of one (1) second or less have little 
effect. 


The effects of a small, sharp turn by the aircraft on 
the position and velocity estimates were found to be fairly 
negligible. On the other hand, the errors in position and 
velocity estimates for an aircraft in a constant turn were 
found to be so large that they would generally preclude the 
use of the present filtering equations (For a turn rate of 
three (3) degrees per second a miss distance of 1100 feet 
would. result, principally from the error in the filtered 
value of ydot.); however, this type of maneuver. should 


never occur during the final precision guidance stage. 


Be. COMMENTS 


A number of comments pertaining to our study should be 
made. First, a number of errors and their sources were not 
Studied, primarily since little could be done at this level 
from a mathematical or a scientific nature to check or 
improve those errors. Such errors include radar alignment 
error, the actual bomb alignment on the aircraft, the 
dampness or dryness of the bomb itself, and the behavior of 


air currents in the atmosphere. 


Proper radar alignment depends upon proper training and 
discipline of the personnel involved in setting up the 
equipment, as well as good surveys prior to setting up the 
radar. These are factors that must be considered in the 
field where the equipment is operated. 


The assumption in this thesis that the tail vanes kept 
the bomb aligned with its axis tangential to the velocity 
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vector relative to the air is felt to be probably the 
greatest source of random error in the system. The bomb 
comes out with angular motion and may even tumble 
momentarily; while this occurs at the time of low speed and 
thinnest air the velocity changes at this time cause large 
integrated effects. Any analysis of this attempting to 
model the aerodynamic forces and moments due to various 
initial configurations would be of doubtful value. The only 
feasible way seems to be numerical experimentation with the 
actual configuration, combined with dropping the bombs and 


observing the paths. 


We also limited ourselves to .one set of standard 
conditions; there is little point to carrying out more 
computations of the same type unless the computations are 
associated with an actual system. Then many of these should 
be carried out over the range of parameters expected to 
occur, such as altitude of plane and of target; speed of the 
plane, and particularly with the best approximation for the 


drag coefficient, as a function of Mach number. 


It must be noted that this research took a general look 
at a sophisticated problem, and in doing so touched only the 
surface in terms of investigating error sources. A  systen 
as expensive and as sophisticated as the AN,yTPQ-27 reguires 
intensive research and testing, as it has received and will 


continue to receive. 


Finally, an important facet of this work was to provide 
a familiarization, or an awareness document for ready 
reference to those individuals who need Some introduction to 
the mathematics used to analyze the system. In that sense 


the intended goal of the research was achieved. 
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